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Abstract. We present the analysis of baryonic and non-baryonic matter distribution in a sample of ten nearby 
clusters (0.03 < z < 0.09) with temperatures between 4.7 and 9.4 keV. These galaxy clusters are studied in detail 
using X-ray data and global physical properties are determined. Correlations between these quantities are analysed 
and compared with the results for distant clusters. We find an interesting correlation between the extent of the 
intra-cluster gas relative to the dark matter distribution. The extent of the gas relative to the extent of the dark 
matter tends to be larger in less massive clusters. This correlation might give us some hints on non-gravitational 
processes in clusters. We do not see evolution in the gas mass fraction out to a redshift of unity. Within rsoo, the 



mean gas mass fraction obtained is (0.16 ± 0.02)/i 



-3/2 
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1. Introduction 

Clusters of galaxies can be regarded in many respects as 
representative for the universe as a whole. As clusters ac- 
cumulate mass from a large volume the baryon fraction 
in clusters is representative for the baryon fraction in the 
universe and hence can be used to determine 17 m when 
comparing the cluster baryon fraction to the upper limit 
of Obaryon from primordial nucleosynthesis. 

The total mass of a cluster can be determined in vari- 
ous independent ways: strong and weak lensing, X-ray ob- 
servations, galaxy velocities and the Sunyaev-Zel'dovich 
effect. From these mass determinations an average baryon 
fraction of 15-20% was inferred (e.g Mohr et al. IT39T31 
Ettori & Fabian [1999 Arnaud & Evrard ITM?1 Grego et 
al. 2001 ) which implies an upper limit of the mean matter 
density O m ^0.3 - 0.4. 

Furthermore, the distribution of the different cluster 
components can be studied when mass profiles are avail- 
able. From X-ray observations the gas density and the 
gas temperature can be determined, which yield not only 
the gas mass profile, but with the assumptions of hy- 
drostatic equilibrium and spherical symmetry also the 
total mass profile of the cluster. Numerical simulations 
showed that these masses are quite reliable for virialised 
clusters (Evrard et al. 11397)1 Schindler [13551 Dolag & 
Schindler 2003 • A comparison of both radial profiles gives 
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information on the relative distribution of dark and bary- 
onic matter. This can be used to learn about the cluster 
formation process: whether all the energy comes only from 
the gravitational collapse, in which case both distributions 
baryonic and dark matter are expected to have the same 
distributions, or whether there are additional physical pro- 
cesses involved. 

Deviations of the Lx — T relation compared to the 
relation expected for self-similar scaling (e.g. Arnaud & 
Evrard 1999) and entropy studies (Ponman et al. 1999) 
suggested that at least for small galaxy systems additional 
(pre-)heating processes (Tozzi |J001) play an important 
role. Numerical simulations were performed to test this 
suggestion (Bialek et al. 120011 Borgani et al. I2001[) finding 
that an entropy floor around 50 - 100 keV/cm 2 is required 
to fit the observational results. 

In this article we determine various cluster properties 
for massive systems and compare the distribution of the 
different components within these clusters. Moreover, we 
use the determined masses to derive relations between the 
masses and other cluster quantities, which give more in- 
formation about cluster formation and evolution. 

We use a sample of nearby clusters with the best 
ROSAT and ASCA data. This sample is complemented 
with four new observations of distant clusters performed 
by Chandra and XMM. 

Throughout this paper we use Ho = 50 km/s/Mpc and 
qo = 0.5. 
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2. The sample 

Since the aim of this paper is the analysis of the total and 
gas mass distribution in nearby galaxy clusters we selected 
clusters in which an accurate total mass determination 
is possible, i.e with relaxed and symmetric morphologies, 
good temperature measurements and well determined sur- 
face brightness profiles. Therefore we obtain a high qual- 
ity sample that consists of ten best clusters in the redshift 
range 0.03 < z < 0.09 observed with the ROSAT PSPC. 
Obviously, the sample is not complete in any sense, there- 
fore no analyses of distribution functions can be made 
with it. But the analyses of correlations between the dif- 
ferent quantities are not affected by the incompleteness 
(see Finoguenov et al. EJ001). Cluster temperatures mea- 
sured by the ASCA satellite are taken from Markevitch et 
al. (HHHHI), Markevitch et al. (HHH, Sarazin et al. {1555^ 
and Bauer & Sarazin ( 2000) resulting in a temperature 
range of 4.7-9.4 keV. In Tabled the different properties 
for each cluster are listed. 

Of all the new observations only the best data are se- 
lected in order to obtain the most reliable X-ray mass 
estimate. In Tabled the published data we use to derive 
masses are listed. 

3. Data analysis 

We use X-ray imaging data retrieved from the ROSAT 
archive 1 to determine the surface brightness profiles of 
the clusters. For each cluster a ROSAT PSPC image was 
reduced using the standard analysis with the EXSAS soft- 
ware. In order to maximize the signal-to-noise ratio, we 
choose the hard energy band (0.5-2.0 keV) corresponding 
to channel numbers 52-201. 

The images were corrected for exposure variations 
and telescope vignetting using exposure maps. The clus- 
ter cores are somewhat blurred in PSPC images by the 
point spread function (PSF) (« 20 - 30 arcsec FWHM for 
ROSAT/PSPC). This effect is more important for com- 
pact clusters. For example for the compact cluster A780, 
Neumann & Arnaud (1999 ) estimated the core radius to 
be overestimated by about 10%. For the other clusters in 
our sample, which have a larger core radius, we estimate 
that the effects of the PSF are much smaller than the 
statistical errors. Therefore no correction for the PSF is 
necessary. 

We generate radial surface brightness profiles in con- 
centric annuli (centred on the emission maximum of the 
cluster) excluding obvious point sources manually. The 
observed profiles are fitted with a /3-model, (Cavaliere & 
Fusco-Fermiano 1976) plus background. 



S(b) = S Q 



-3/3+i 



B. 



(1) 



The parameters r c (core radius), /3, Sq (central sur- 
face brightness) and the background B are obtained from 

1 http:/ / www.xray.mpe. mpg.de/rosat/archive 



a least-squares fit to the X-ray profile. The slope f3 and 
the core radius r c are not independent parameters, with 
/3 increasing when r c increases. They are found to range 
between 0.6 and 0.8 and between 130 and 290 kpc, respec- 
tively. 

However, the overall /3-model fit is a poor description 
of the central region of some clusters where excess emis- 
sion is observed. Indeed we find in most of the cases very 
large \ 2 values when fitting the entire cluster emission. We 
reduce the x 2 values by excluding the central bins from 
the fit. The best fit /3 -model was determined excluding 
the data within the cooling radius taken from Peres et 
al. C2S3, Allen & Fabian flggj , White et al. itTMTIl . 
For the clusters in our sample with excess central emis- 
sion, the exclusion of the central part of the profile yields 
larger f3 values and core radii values compared to the over- 
all 0- model fit. Excluding the central excess we underes- 
timate the gas density at the centre (about 12% for the 
central 3' in the case of the cluster A2199). However we 
estimate that the central gas mass contributes only with a 
few percent in the gas mass at larger radius (about 3% for 
cluster A2199 at the radius of r^oo). Therefore for the gas 
mass determination, this underestimate at small radii is 
negligible when we integrate the gas density out to large 
radii. 

In Table |21 the resulting parameters are shown. The 
reported errors are 90% confidence level. 

4. Mass determination 

Once we deproject the surface brightness to three- 
dimensional density with the f3 model, together with the 
assumptions of hydrostatic equilibrium and spherical sym- 
metry the integrated mass within radius r can be deter- 
mined as 



[im v G 



din; 



d\ni 



(2) 



with p and T being the density and the temperature of 
the intra-cluster gas, respectively, k, fi, m p , and G are 
the Boltzmann constant, the molecular weight, the proton 
mass, and the gravitational constant. 

4.1. Isothermal analysis 

As well as the density profiles we need temperature pro- 
files to determine the total cluster mass. We follow two dif- 
ferent approaches for the temperature profile. In our first 
study we assume isothermality. In Sect. 14.21 we include the 
temperature gradients derived from Markevitch's ASCA 
analysis to see how the total cluster mass is affected. In 
the isothermal approach we neglect the term associated 
with the temperature gradient in Eq. J2J). In this case the 
total mass profile is: 



M(< r) = 



3/cT t 



gas' 



1 



r), 



G/im p r 2 + ri 
where the mean atomic weight p is assumed to be 0.61 



(3) 
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Table 1. In Cols. 1 and 2 the clusters and their redshift are shown. In Col. 3 the hydrogen column density from Dickey 
& Lockman (1990) is listed. The temperatures in Col. 4 are ASCA emission- weighted temperatures excluding cooling 
flows and other contaminating components. Col. 5 lists the ratio of the bolometric flux to the ROSAT PSPC countrate in 
the energy range 0.5-2.0 keV using the programme PIMMS (http:/ /www. heasarc.gsfc.nasa.gov/Tools/w3pimms. html). 
ASCA temperatures are taken from references: 1— Markevitch et al. (1998); 2— Markevitch et al. (1999); 3— Sarazin et 
al. (1998); 4-Bauer & Sarazin (2000). 
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A644 
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4 
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Table 2. Published data from XMM and Chandra observations used to derived the total and gas mass. In Cols. 1 and 
2 the clusters and their redshift are listed. The gas temperature is shown in Col. 3. The published fit parameters of the 
(3 model are listed in Cols. 4, 5 and 6: the slope /3, the core radius r c and the central electron density n c o. The quoted 
errors are 90% confidence level for all quantities except for the parameters in cluster RX-J0849+4452 and the value 
of the temperature in A1835 which are 68% confidence level. References: (1) Standford et al. ( 2001); (2) Schindler et 
al. ijMFl) : (3) Arnaud et al. lj2T)I72^ : (4) Majerovitch et al. plH^ . 



Cluster 


z 


(keV) 


P 


(kpc) 


(lO^cm -3 ) 


Ref. 


RX-J0849+4452 

RBS797 
RX-J1120.1+4318 
A1835 


1.26 
0.354 
0.6 
0.25 


5.8it 8 7 

7 7+1.2 
' • ' -1.0 

5.3 ±0.5 
7.6 ±0.4 


0.61 ±0.12 
0.63 ±0.01 

n 70+O.O6 

0.704 ± 0.005 


100.2 ± 30.7 
49 ±4 
20lt£ 

202.3 ± 7.1 


1.42 
8.86 
0.81 
1.47 


Chandra (1) 
Chandra (2) 
XMM (3) 
XMM (4) 



The total mass thus depends linearly on both (3 and 
Tgas- In the case of clusters with central X-ray excess, we 
use emission-weighted gas temperature obtained by ex- 
cluding the central part of the cluster. Typical mass pro- 
files are shown in Fig. ^ for the cluster Abell 2199. 

After having determined the gravitational mass pro- 
files for the clusters, it is important to fix the radius 
within which the cluster masses can be calculated. As 
the mass of a cluster increases with radius, masses can 
only be compared when derived within equivalent vol- 
umes. Simulations by Evrard et al. ( 1996) showed that the 
assumption of hydrostatic equilibrium is generally valid 
within at least radius r5oo, where the mean gravitational 
mass density is equal to 500 times the critical density 
Pc(z) = 3i?Q (1 ± z) 3 /8nG. The resulting cluster masses 
and gas masses within r5oo and 0.5 x rsoo for all the clus- 
ters are listed in Table 0] 

The errors on the total mass can be estimated as fol- 
lows. The total cluster mass is affected by the errors in 



the parameters of the /3-model. This error is estimated 
to be about 5%. The uncertainties in the total mass esti- 
mate are much larger due to the uncertainty in temper- 
ature estimates and possible temperature gradients. For 
larger radii we assume an error of 10% in the total cluster 
mass, caused by the existence of a temperature gradient. 
There are also additional uncertainties coming from devi- 
ations from spherical symmetry (Piffaretti et al. 2003), 
deviations from hydrostatic equilibrium and projection 
effects (together about 15% - 30%, Evrard et al. ITM61 
Schindler I1996f) . but these are hard to quantify for each 
cluster individually. As only well relaxed clusters where 
chosen, these errors should be relatively small (< 15%) in 
the clusters of this sample. 

We compute the error in the total cluster mass as the 
convolution of the error coming from the uncertainty in 
the fit parameters, the error in the temperature and the 
10% error coming from the assumption of non isothermal- 
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Table 3. X-ray quantities as measured from ROSAT/PSPC and ASCA observations. The clusters are listed in Col. 1. 
Col. 2 shows the emission- weighted gas temperature. In Cols. 3, 4, 5 and 6 the fit parameters of the f3 model are shown: 
the slope /3, the core radius r c , the central surface brightness 5*0 and the background B in the energy band 0.5 - 2.0 
keV. So is in units of 10~ 2 ROSAT/PSPC counts/s/arcmin 2 . B is in units of 10" 4 ROSAT/PSPC counts/s/arcmin 2 . 
In Col. 7 rf denotes the radius range fitted. 
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8.0±0.4 
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5.4 


2.5 


0.11 - 
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A478 


8.4±;-s 


0.69 ±0.02 


201 ± 20 


11.7 


1.5 


0.22 - 


6.1 


A644 


8-6±°iS 


0.72 ±0.02 


239 ± 17 


6.9 


2.3 
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5.4 


A2029 
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0.68 ± 0.02 


244 ± 24 


9.7 
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Table 4. Results of the isothermal analysis: total mass, gas mass and gas mass fraction for the cluster sample. The 
first column gives the cluster name. Col. 2 denotes the radius t^qq which comprises an overdensity of 500 over the 
critical density. Cols. 3, 4 and 5 list the total mass, the gas mass and the gas mass fraction within rsoo, respectively. 
In Cols. 6, 7 and 8 the same quantities are listed for a radius 0.5 x r^oo. For all the quantities the reported errors are 
90% confidence level. 



Cluster rsoo 


MtotfVsoo) 


M gas (r 5 oo) /gas( r 50o) 




M gas (IfH) /gas(^f) 


(Mpc) 


(10 14 M Q ) 


(10 14 M Q ) 


(1O 14 M ) 


(10 14 M Q ) 



A496 


1 


.42 ± 0.03 


4.6 ±0.5 


0.69 


A2199 


1 


.45 ± 0.03 


4.8 ±0.6 


0.68 


A3112 


1 


44+0-09 

'^-0.14 


5-3±S;i 


0.94 


A1651 


1 


.55 ±0.05 


6.9 ±0.9 


1.24 


A3571 


1 


7 q + 0.03 

• '°-0.02 


8.5 ± 1.0 


1.50 


A1795 


1 


on+0.12 
• 8 ^-0.13 


10.4 ± 1.8 


1.49 


A401 


1 


7O+0.04 
•' 8 -0.05 


10.2 ± 1.2 
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A478 
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•""-■-o.ie 


11.5±J;5 


1.83 


A644 


1 


.94 ± 0.07 


13.0±if 


1.52 


A2029 


1 


nr+0.06 
Md -0.05 


13.4 ± 1.7 


2.16 



0.150 ±0.019 2.2 ±0.3 0.27 0.12 ±0.03 

0.142 ±0.018 2.4 ±0.3 0.28 0.12 ±0.03 

0-181^3 2.5 + _° .t 0.37 0.15±g;82 

0.18 ±0.02 3.2 ±0.4 0.50 0.16 ±0.04 

0.18 ±0.02 4.0 ± 0.5 0.60 0.15 ±0.04 

0.14 ±0.02 5.0 ± 0.9 0.64 0.13 ±0.05 

0.19 ±0.02 4.7 ± 0.6 0.75 0.16 ±0.04 

0.W±° o ° l 5-5±S1 0.80 0.U±° o ° i 

0.117 ±0.017 6-2±° 9 s 0.68 0.11 ±0.03 

0.16 ±0.02 6.4 ± 0.8 0.90 0.14 ±0.04 



ity However we keep in mind that the true uncertainty of 
the total mass is probably greater than this value. 

In the energy range considered, the gas mass estimate 
is almost independent on the temperature measurement. 
We estimate the error in the gas mass to about 5% due to 
the uncertainty in the fit parameters. 

4.2. Temperature gradients 

The isothermal assumption may give poor estimates of the 
total mass if strong temperature gradients are present. 
Observations with ASCA suggest that the temperature 
does decrease with radius. To estimate the effect of these 
temperature gradients on cluster masses, we have cal- 
culated Mtot using the temperature profiles derived by 



Markevitch et al. 119981 Markevitch et al. 119991 Sarazin et 
al. 119981 and Bauer & Sarazin ESfjfjQI In some clusters, the 
temperature outside the cooling core can be well approx- 
imated by a polytropic function 



Tgas °C Pgas 



(4) 



where 7 is the polytropic index. In this case the total mass 
enclosed in a sphere of radius r (Eq. 0) is: 



M(< r) = 



3/3 7 ^ 3 kT gas (r) 
Gm v [i (r 2 + r 2 ) 



(5) 



Here T gas is the true temperature, rather than a projection 



on to the plane of the sky. Markevitch et al. l|1999f) showed 
that as long as the temperature is proportional to a power 
of the density, and the density follows a /3-model, the true 
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Fig. 1. Total mass profile assuming isothermality (solid 
line) and gas mass profile (dotted line) for the cluster 
A2199. The dashed line represents the total mass profile 
derived with a temperature gradient. r 500 is the radius 
which comprises an overdensity of 500 over the critical 
density, represents the maximum radius out to which 
the temperature profile is calculated. A typical error bar 
for the isothermal total cluster mass is shown at r^oo. 

temperature differs from the projected temperature only 
by a constant factor, given by: 

Jptoj = r[f/?(i + 7 )-|]r(3/?) 

r truc r[|/3(i + 7 )]r(3/3-i)- (> 

In our analysis, the projected temperature profiles are 
parametrized using a linear function of the form 

T gas (r) = T(0) - or. (7) 

The parameters T(0) and a are determined by fitting a 
straight line to the projected temperature profiles men- 
tioned before, excluding the central cluster region in the 
case of clusters with central X-ray excess. It is not useful 
to fit a more complex function because the temperatures 
for consecutive annuli are determined with low accuracy. 

We do not deprojected the temperatures for the mass 
analysis with temperature gradient. Although this intro- 
duces some additional uncertainty in the total mass we 
show below that the results using true temperatures or 
projected temperatures are very similar. 

For the cluster A2199 we compare the total mass 
derived using the linear gradient Eq. Q and the mass 
calculated with the polytrope Eq. (0J by Markevitch et 



Fig. 2. Gas mass fraction assuming isothermality (solid 
line) for cluster A2199. Dashed line represents the gas 
mass fraction derived with a gradient of temperature. A 
typical error bar for the isothermal gas mass fraction is 
shown at rsoo- 

al. ( 1999). The estimated values for the total cluster mass 
agree within the errors. There is a difference < 10% at 
large radii (r > IMpc) and < 15% at the small radius 
(r = Q.2Mpc). For example Markevitch et al. I|1999[) using 
a polytrope equation with 7 = 1.17 find a total cluster 
mass of (0.65 ±0.11) x 10 14 Af o , (2.9 ±0.3) x 10 14 M Q and 
(3.6 ± 0.5) x 10 14 M o at a radius of 0.2 Mpc, 1 Mpc and 
^500 = 1.3 Mpc respectively, for cluster A2199. We derive 
the total cluster mass for this cluster using the linear gra- 
dient T(r) = — 2.D' + 5.6 where r is units of Mpc and tem- 
perature is in keV. The masses obtained in this case are 
quite similar, (0.56±0.14) x 10 14 M o , (3.1±0.8) x 10 14 M Q 
and (3.9 ± 1.0) x 10 14 M o at a radius of 0.2 Mpc, 1 Mpc 
and 1.3 Mpc, respectively. Therefore we conclude that for 
our purpose of comparison, it is sufficient to use a linear 
gradient of projected temperatures. 

In the following we will apply only a linear temperature 
gradient. In this case we estimate the errors in the total 
cluster mass and gas mass as follows. Due to the errors 
in the temperature, the errors in the linear fit are large. 
We estimate the errors in the total cluster mass using the 
different possible temperature gradients given by the lower 
and upper values of T(0) and a (see Eq. Q). When the 
different possible temperature gradients are used, the radii 
7Y5oo and 7Y500 /2 change significantly and hence the total 
and gas mass enclosed by them. For the gas mass an error 
of ~ 10% is estimated. 
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Fig. 3. Comparison of masses obtained with the isother- 
mal and temperature profile analyses at radius 7"5oo- The 
smaller symbols represent the temperature profile analy- 
sis. 

In Fig. ^ t ne different mass profiles calculated using 
the isothermal model, and a linear gradient of tempera- 
ture, for cluster A2199 can be seen. With the assumption 
of constant temperature (solid line) the total mass is over- 
estimated (~ 10%) at large radii and underestimated at 
small radii, compared to the temperature gradient anal- 
ysis (dashed line). The overestimate of the total mass is 
only significant at radii larger than n, with ri being the 
radius out to which the temperature was measured. This 
trend is observed for all the clusters in our sample. Fig. ^ 
also shows the gas mass profile (dotted line) in the radius 
range that has been used to fit the data with the /3-model. 

Fig. |21 shows the gas mass fraction, defined as the 
ratio between the gas mass M sas and the total mass 
M tot . The difference in the gas mass fraction between the 
isothermal and non-isothermal cluster studies is shown 
as well (solid and dashed lines, respectively). At small 
radii both analyses provide approximately parallel profiles 
with smaller values obtained in the temperature gradient 
case. At larger radii the temperature profile analysis yields 
steeper profiles and higher values (~ 10%) compared to 
the isothermal one. The gas mass fraction derived using a 
gradient of temperature is not reliable beyond the radius 
tt where the gradient is determined and thus it is not 
plotted. 

The masses calculated with both analyses (isothermal 
and temperature profile) are compared in Fig. for each 
cluster in the sample at radius rsoo- The symbols corre- 



Fig. 4. Gas mass fraction profiles derived for the nearby 
cluster sample. Profiles are plotted from the minimum ra- 
dius fitted in the (3 model. 

spond to the cluster name listed in Fig. The smaller 
symbols represent the values for the temperature profile 
analysis. Although the gas mass profile is not influenced 
by the change in the cluster temperature, the gas mass 
at radius rsoo is different in the two analyses due to the 
difference in rsoo which depends on the temperature. 

5. Results and discussion 

In the following all the results presented are refered to the 
analysis of the cluster sample assuming isothermality. 

5.1. Gas mass fraction 

The gas mass fraction defined as f gas = M gas /M to t is cal- 
culated for each cluster in the sample. The errors asso- 
ciated with the gas mass fractions are calculated by the 
convolution of errors in the total cluster mass and the gas 
mass. We find that inside each cluster the gas mass frac- 
tion increases with radius (see Fig. implying that the 
gas distribution is more extended than dark matter. 

A low fi is required to reconcile the high gas mass 
fraction of < / gas > r500 = 0.16 ± 0.02, whith the baryon 
fraction predicted by primordial nucleosynthesis. 

To test whether there is any evolution of the gas 
mass fraction we plot this quantity versus redshift 
(Fig. |5j) including the analysis from Schindler (1999) 
for distant clusters (0.3<z<1.0). In this figure, as 
well as in all other following figures, each cluster 
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Fig. 5. Gas mass fraction versus redshift. No clear evolu- 
tion of the gas mass fraction can be seen out to a redshift 
of unity. See Fig. for the different cluster symbols. 
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Fig. 6. Symbols used for the various clusters in the figures. 

is plotted with a different symbol (see Fig. 0). We 
include in the comparison another distant clusters: 
RX-J0849+4452 (z=1.26), RBS797 (z=0.354), A1835 
(z=0.25) and RX-J1120.1+4318 (z=0.6) where we calcu- 
late the total and gas mass using the published parameters 
from the Chandra and XMM data analysis (sec Table [2J . 

In the gas mass fraction we see no clear trend with red- 
shift. The mean value in our sample is < / gas >= 0.16 ± 
0.02 which is in agreement with other nearby samples, 
Mohr et al. (1999): / gas =0.21, Ettori & Fabian (1999): 



Fig. 7. Gas mass versus total mass. In solid line is shown 
the best fit for the nearby sample. For comparison, in dot- 
ted line is plotted the best fit obtained for the distant 
sample by Schindler (1999) 

/gas=0.17 and also in agreement with the value for distant 
sample by Schindler (1999), / gas =0.18. Therefore we con- 
clude that we do not see evolution in the gas mass fraction 
out to a redshift of unity. In contrast to this result Ettori 
& Fabian (1999) found indications for a decrease of / gas 
with increasing redshift in a nearby sample. Matsumoto 
et al. (20O found no clear evidence of evolution of / gas 
for the clusters at z < 1.0. 

In Fig. [7] we compare the gas mass with the total clus- 
ter mass. A trend of an increasing gas mass with total mass 
is visible in our sample. This trend was also found for dis- 
tant clusters by Schindler l|1999|l (dotted line in Fig. 0). 
A linear fit regression yields 

M gas (r 500 ) = 0.13M tot (r 50 o) (1 - 09±al2) . (8) 

M gas (r 5QQ ) and M tot (r 50 o) are in units of 10 14 M Q . This 
trend is in agreement with the analysis by Arnaud & 
Evrard <HH2|. 

Because the exponent in Eq. (JSJ) is close to unity, we 
find no clear dependence of the gas mass fraction on the 
total mass in our sample (see Fig. [8] where the nearby and 
distant sample are shown). 

Several authors have also related measured values of 
/gas and M to t or T gas . Up to now the results are discordant 
concerning the form of this relation. For example, David 
et al. ( 1995) found indications for an increase of / gas with 
increasing T gas . Allen & Fabian l|1997f) found in a sample 
of X-ray luminous clusters indications for a decrease of 
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Fig. 8. Gas mass fraction versus total cluster mass for the 
nearby and distant sample. The gas mass fraction shows 
no clear dependence on the total mass. 



/g as with increasing T gas . Ettori & Fabian (^999) found no 
dependence of / gas on M to t for high luminosity clusters. 

As mentioned before the gas mass fraction is not con- 
stant with radius. In order to plot this increase against 
other cluster properties, we compare the gas mass frac- 
tion at radius r 50 o with the gas mass fraction at 0.5 x r 500 
in each cluster. The mean gas mass fraction at 0.5 x r^oo 
is 0.138 ± 0.016, i.e. smaller than the mean of 0.16 ± 0.02 
at rsoo- The ratio of these fractions is a measure for the 
extent of the gas distribution relative to the dark matter 



extent: E 



/gas (^500 ) 



Wrsoo/2)- 

The error in the relative gas extent E is not easy to de- 
termine. We estimate its uncertainty by testing how much 
this value changes when: 

• we consider the uncertainty in the gas temperature 

• the fit parameters change 

Since the relative gas extent is calculated as the ratio 
of gas mass fractions at different radii, the errors coming 
from the uncertainty in the temperature are cancelled out. 
The relative gas extent E changes significantly when the 
fit parameters (/? and r c ) are varied. This error coming 
from the uncertainties in the fit parameters is included in 
the following graphs. 

For all the clusters in our sample the relative gas ex- 
tent E is larger than 1 (see Fig. |5J). This means that in 
general the gas distribution is more extended than the 
dark matter, which is in agreement with other results for 
nearby cluster samples by David et al. (|1995[1 . Jones & 



Fig. 9. Ratio of gas mass fraction at rsoo arid rsoo/2 as 
a measure for the relative extent of the gas distribution 
versus total cluster mass at r^oo- For clarity, clusters with 
the largest error bars from Schindler ( 1999) are not shown. 



Forman (1999), Ettori & Fabian (199SJ) and for distant 
samples by Schindler l(T555)l and Tsuru et al. 1(1555)1 . The 
distant sample by Schindler ( 1999) and the other Chandra 
and XMM observations are also shown in Fig. |5J for com- 
parison. 

In the nearby sample (4.7 < T < 9.4 keV) this relative 
gas extent E shows a mild dependence on the total cluster 
mass (see Fig.[5| similar to the result by Schindler JT999). 
Clusters with larger masses tend to have smaller relative 
gas extents (similar dependence confirmed by Reiprich & 
Bohringer 119991 although they used different radii). 

The differences in the gas and dark matter distribu- 
tion cannot be explained by purely gravitational energy 
coming from the collapse of the cluster. Additional energy 
input is necessary to explain it, e.g. from supernova-driven 
galactic winds. It might be that this additional energy af- 
fects low mass clusters more than massive clusters, so that 
a massive cluster can maintain a ratio E ss 1, while in the 
smaller clusters the gas is more extended. 

Entropy studies by Ponman et al. (1999) of cool clus- 
ters (T < 4keV) observed with ROSAT and GINGA sug- 
gested that for these systems, (pre-)heating processes play 
an important role in cluster formation. 

Equation @ shows that for a given radius M to t oc 
T x (3. If there is a temperature rise due to preheating 
processes, a cluster with a certain M to t should have shal- 
lower gas density profile or more extended gas distribution 
{(3 value smaller). Then if that heating is more effective 
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Fig. 10. /3-temperature relation. (3 is the parameter ob- 
tained when fitting with a (3— model the X-ray surface 
brightness profile. There is a trend to find larger (3 values 
for larger temperatures. 



in cooler cluster, i.e. less massive clusters, one should ex- 
pect an anticorrelation between total mass and gas extent 
and lower values of [3 in cooler systems. This seems to be 
consistent with observations (e.g. Mohr & Evrard 119971 
Vikhlinin et al. [1999, Ponman et al. UMl Nevalainen 
et al. I2UUU|) . Our sample also exhibits this behavior, see 
Fig.Hni 

5.2. Mass-temperature relation 

Accurate measurements of the cluster total mass are pos- 
sible only for a limited number of clusters. For this reason 
there arc currently not enough data available for a direct 
derivation of the mass function, which is crucial for the 
determinations of the cosmological parameters using the 
cluster abundances at different redshifts. A more practical 
way of determining the mass function is to observe the dis- 
tribution of readily available average cluster gas temper- 
atures and to convert these to masses, taking advantage 
of the tight correlation between mass and temperature 
predicted by hydrodynamic cluster formation simulations 
(e.g. Evrard et al. I1996|l . Therefore, a well-established 
Mtot — T relation can be used as a powerful cluster mass 
estimator. Furthermore the M tot — T relation is also in- 
teresting in itself, because deviations from the predicted 
self-similar scaling of Al oc T 3 / 2 would indicate that more 
physical processes are at play than gravity alone. 



Assuming self-similarity and a velocity dispersion pro- 
portional to the X-ray temperature, the virial theo- 
rem provides a relation between total mass, radius and 
X-ray temperature: M to t('"5oo)/''500 oc T. Equivalently, 
f5oo can be expressed by the definition of the overden- 
sity r 50 o oc M to t(r 5 oo) 1/3 /(l + z) yielding the relation 
M to t(f*5oo) oc (T/(l + z)) 3 / 2 . We see this correlation in 
our data (see Fig.lllb). A fit taking into account the errors 
in the temperature and the error in the mass yields the 

/ rp n(1-7±0.2) 

following relation Mt t(^50o) = 0-36 1 jqr^ J shown 

in Fig. Illb as solid line. Mtot^soo) is in units of 10 14 M Q 
and T in units of keV. The slope is greater than the virial 
value of 3/2 but consistent within the errors. If we do the 
same fit using the masses calculated with the temperature 
profiles the slope found is 1.4 ± 0.2. 

Other observational analyses, with only high- 
temperature (kT > 4 — 5keV) clusters, have achieved 
results consistent with the theoretical prediction (Hjorth 
et al. HTKM Neumann & Arnaud HUH . 

However many studies (e.g. Ponman et al. 1999 
Horner et al. 1999, Nevalainen et al. 2000, Finoguenov et 
al. l2U01[> have shown that the influence of energy feed-back 
into ICM in low-temperature clusters/groups can become 
more significant than that in high-temperature systems. 
Consequently, the self-similarity may break at the low- 
temperature end. If supernovae release a similar amount 
of energy per unit gas mass in hot and cool clusters, the 
coolest clusters would be affected more significantly and 
exhibit a stronger shift to higher temperatures in the M- 
T diagram than the hotter clusters. This will steepen the 
M-T relation. 

Furthermore, we test the relation between the gas 
mass and gas mass fraction and the temperature. 
As expected from the gas mass— total mass relation, 
there is also a relation between gas mass and tem- 
perature (see Fig. II lb) . A linear regression fit yields 
M gas (r 50 o) = O.OITt 1 - 80 "- 16 ) with M gas (r 50 o) in 'units 
of 1O 14 M0 and T in keV. This correlation is in agree- 
ment with other results in nearby clusters (Reiprich & 
BohringerHnHHI Jones & Forman [1.999). For comparison, 
Reiprich & Bohringer (1999) find an exponent of 2.08 and 
Schindler (1999) a larger exponent 4.1±1.5 for distant 
clusters. 

As expected from the non-correlation of the gas mass 
fraction with the total mass, we find also no correlation 
between the gas mass fraction and the temperature (see 
Fig. II lb ) . This result is in good agreement with Mohr et 
al. (1999). They find a mild dependence comparing low 
temperature clusters (T < 5keV) with high temperature 
clusters (T > 5keV). For the high temperature clusters 
alone, in which category all our clusters fall, they find no 
dependence. 

We find an interesting correlation between the relative 
gas extent and the temperature (see Fig.lllH). which is of 
course related to the dependence of the relative gas extent 
on the total mass, shown above. The relative gas extent 
tends to be slightly larger in lower temperature clusters. 
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Fig. 11. Various quantities versus temperature, (a) Total mass, (b) Gas mass. Solid lines show the best fit for our 
sample. The dashed line is not a fit but the curve expected from virial considerations with an arbitrary normalisation, 
(c) Gas mass fraction, (d) Ratio of gas mass fractions at rsoo and rsoo/2. 



6. Summary 

We have analysed a sample of ten nearby clusters of galax- 
ies using the X-ray data provided by the ROSAT and 
ASCA satellites. For this sample physical quantities like 
gas mass, total mass, gas mass fraction and relative gas 
extent have been derived. Correlations between the above 
quantities have been studied and our findings are 

• The gas mass fraction increases with radius for all our 
clusters implying that gas is more extended than dark 
matter, confirming previous results (David et al.[l995 



Ettori & Fabian HMfl Jones & Forman[l999 ) . This be- 
haviour is more pronounced when temperature profiles 
are taking into account in the mass analysis. 

• Within rsoo, the mean gas mass fraction obtained is 
(0.16±0.02)/i~ 3/2 . 

We see no trend in the gas mass fraction with redshift. 

• The gas extent relative to the dark matter distri- 
bution shows a mild dependence on the total mass 
and gas temperature. Clusters with larger masses have 
smaller relative gas extents, as we would expect if non- 
gravitational processes are important in cluster forma- 
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tion. Hints for this kind of behaviour have been found 
previously in distant clusters. These new results con- 
firm this trend. 
• Studying the mass-temperature relation we find a slope 
slightly steeper compared with the theoretical value of 
3/2 although consistent within the errors. The self- 
similar slope 3/2 is found when using temperature 
gradient analysis. Other observational analyses, with 
only high-temperature clusters, have achieved results 
consistent with the theoretical prediction (Hjorth et 
al. irgggl Neumann & Arnaud \T$M Nevalainen et 

ai.mp. 
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